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We show that the phase sensitivity A9 of a Mach-Zehnder interferometer fed by a coherent 
state in one input port and squeezed-vacuum in the other one is i) independent from the true 
value of the phase shift and ii) can reach the Heisenberg limit AO ~ 1/Nt, where Nt is the 
average number of particles of the input states. We also show that the Cramer-Rao lower bound, 
A6 oc 1/V|a| 2 e 2r +sinh 2 r, can be saturated for arbitrary values of the squeezing parameter r and 
the amplitude of the coherent mode \a\ by a Bayesian phase inference protocol. 
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Introduction. The goal of quantum interferometry is to 
estimate phases beyond the shot-noise ("standard quan- 
tum") limit. The quest requires proper non-classical 
states, as was first shown by Caves in 1981 [l|, who 
considered a Mach-Zehnder (MZ) fed by coherent <X> 
squeezed-vacuum light. This benchmark generated a 
large body of theoretical [1, H, 0, and experimental 
@, 0, 0] studies, including the demonstration of sub shot- 
noise sensitivity Q using parametric down-conversion 
in a cavity as a source of squeezed vacuum [8]. The 
scheme proposed by Caves is sketched in Fig.([T]). One 
of the inputs of the linear loss-less MZ is the coher- 
ent state \a) a = X)m=o ^ m l m ) a ' w * tn a = e^M an d 
=y= — . The second input state is the squeezed- 



with the average number of photons injected in the MZ 
n ~ \a\ 2 . 

In this Letter we show that the choice of the average 
relative number of photons as phase estimator is not opti- 
mal. Quantum fluctuations also contains information on 
the true value of the phase shift, which can be retrieved 
by taking in account the higher moments of the measured 
number of particles at the output ports. We will show 
that the ultimate phase sensitivity of a Mach-Zehnder 
fed by coherent £g> squeezed-vacuum light is 
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m)f,, with £ = re 3 and 
H m (0) [9j, H m (x) being the Hermite 
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polynomials. Following the current literature, based on 
the original works of the 80's, the phase estimate of this 
system is retrieved from the measurement of the relative 
number of particles at the output ports M out = N c — N^. 
Fluctuations on the results obtained in p independent 
measurements propagate to the estimated value of the 
phase shift 9 |10j | , w hich can be eventually determined 
with uncertainty 111 ]: 
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According to Eq.fTJ, we can appreciate an increase of 
phase sensitivity with respect to the shot-noise only when 
the true value of the phase shift is sufficiently close to 
8 = tt/2 P,i (dark fringe) , where (M out ) = (N c — N d ) = 
0. On the other hand, (M out ) depends weakly on the 
phase shift when 9 w 0, it and the error propagation 
formula Eq.|T]) predicts large phase fluctuations around 
these points. Asymptotically in the amplitude of the co- 
herent state, \a\ 2 3> sinh 2 r, and for a fixed squeezing 
parameter r, Eq.JT]) predicts a sub shot-noise sensitivity 
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The phase sensitivity Eq.([3]) is i) independent from the 
true value of the phase shift over the whole interval < 
< 7T and ii) it reaches, at the optimal point \a\ 2 = 
sinh 2 r, the Heisenberg limit: 
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asymptotically in the average number of photons n — 
\a\ 2 + sinh 2 r and with a number of independent mea- 
surements p > 30. 

In the following, we will first analytically calculate 
the Cramer- Rao lower bound (CRLB), Eq.([3]), and then 
demonstrate that it is saturated by a Bayesian phase 
inference approach. A proof of principle of Eq.Q can 
be obtained within current technology, at least in the 
limit of small n: high-efficiency number-resolving photo- 
detectors have been recently applied to interferometry 
13l . lij and high squeezing has been obtained with para- 
metric down-conversion hj. Out results can be relevant, 
for instance, to improve the efficiency of the large scale in- 
terferometers dedicated to the detection of gravitational 
waves [l6j . which would not require phase-stabilization 
techniques to lock at the optimal point 9 = tt/2 Q and 
which can significantly increase their sensitivity. 

The Cramer-Rao lower bound. The output state 
of a loss-less Mach-Zehnder interferometer is given by 
\ipout) = e~ ie J y \ipi n ) [13], where, in our case, \tpi n ) = 
\&}a\Qb- The conditional probability to measure N c and 
Nd particles at the output ports, given an unknown phase 
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FIG. 1: (Color online) Schematic representation of the Mach- 
Zehnder interferometer. The input modes a, b are a coherent 
and a squeezed- vacuum field, respectively. 



shift 9, is 



P(N c ,N d \9) = 



N 
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where /Lt = (N c - N d )/2 
matrix elements. The 



and d^„(cos9) are 
Fisher information, 

dP(N c ,N d \e)\2 



rotation 
F(9) = 

l^N a =0 2^N d =Q P(N c ,N d \8) \ de J , turns oui 10 oe 
independent from the true value of the phase shift 9, 
see Fig.©, and an analytical calculation gives F{9) = 
\a\ 2 e 2r + sinh 2 r. According to Cramer and Rao, the 
phase sensitivity of an unbiased estimator is bounded by 
A9 = . 1 , which, after replacing the previous expres- 

y'pF(S) 

sion for the Fisher information, gives Eq.©. There are 
interesting limit regimes recovered by this equation: i) 
When r = or a — we get the (^-independent) shot- 
noise limit A9 — 1/ y/pn. The phase independence of the 
case r = has been studied and experimentally demon- 
strated in [l3j |. ii) When sinh 2 r <C \a\ 2 we obtain the 
sub shot-noise limit discussed by Caves, A9 = e~ r / y/jm 
Eq.© with, again, the important difference that, here, 
the phase sensitivity is independent of the phase shift for 
^ 9 ^ 7r. Notice that, in the limit of very high squeez- 
ing, sinh 2 r > |a| 2 , Eq.© predicts A9 = 1/y/pn (still 
at 6 = 7r/2), while Eq.© gives a sub shot-noise scaling 
A9 = l/(v^ v / 4japTl). 

The most important regime predicted by Eq.© is ob- 
tained when \a\ 2 ~ sinh 2 r — n/2 (i.e., with half of the 
input intensity provided by the coherent state and half 
by the squeezed light). This gives A9 = when 
n,p^S> 1. It is interesting to notice that, for these optimal 
values of the parameters a and r, the error propagation 
formula Eq.© predicts a divergence. In figure ©a) we 
compare, as a function of r and for 9 = 7r/2, the quantity 
y/rvp A9 calculated with Eq.© (dotted line) compared 
with Eq.© (solid line). 

Why does the error propagation formula Eq.© pro- 
vides such a poor phase sensitivity with respect to the 
CRLB? The answer is that an estimation of the phase 
shift based only on the measurement of the average rela- 
tive number of particles does not exploit all the available 
information contained in the detection of N c and N d . It 
forgets about the information contained in the fluctua- 



FIG. 2: (Color online) Comparison between Eq.(JTJ (dashed 
red line) and Eq.© (solid blue line). Circles are the results of 
the Bayesian analysis. A) Phase sensitivity A8^/np as func- 
tion of the squeezing parameter r, for 8 = n/2 and \a\ 2 = 10. 
Notice that Eq.JTJ diverges at sinh 2 r = \a\ 2 (dotted verti- 
cal line). For r » 1, Eq.© gives A6^/np 
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(dotted horizontal line). B) AO^/n/p in the limit p — *■ oo as a 
function of the true value of the phase shift. Here \a\ 2 — 10 
and r = 1. 



tions of the total number of particles and (since the rel- 
evant probability distributions are not Gaussians) in the 
higher moments [20|. In figure ©b) we plot the phase 
sensitivity s/np A9 as a function of the true value of the 
phase shift. The dashed line is the result Eq.© and the 
solid line is Eq.©. 

Bayesian analysis. Is it possible to saturate the CRLB 
and demonstrate a phase sensitivity at the Heisenberg 
limit A9 ~ I/Nt, being Nt the average number of 
particles burnt during the estimation process [2l|? A 
possibility, of course, is to consider the Maximum Like- 
lihood estimator which, according to the Fisher theo- 
rem, saturates the CRLB asymptotically in the num- 
ber of measurements p. In the following, however, we 
consider a Bayesian protocol [19( showing that it also 
saturates the CRLB. To simulate a phase estimation ex- 
periment, we i) randomly choose p values Nc , at 
the output ports distributed according to P(N C , Nd\0) 
with an unknown 9; ii) invert, by applying the Bayes 
theorem, the distribution Eq.© and associate to the 
measured values N^ , N^ the probability distribution 
P(<f>\N^\NP;...;N^,N^) ~ U.U , N^); 

iii) calculate the phase sensitivity as 68% confidence 
around the maximum of the phase distribution. In 
Fig.©a,b), the circles, obtained with the Bayesian prob- 
abilities asymptotically in the number of independent 
measurements p, coincide with the analytical expression 
of the CRLB, Eq.©. 

Yet, in order to demonstrate the possibility to reach 
the Heisenberg limit, A9 ~ 1/iVr, we have to carefully 
analyze the role of p [13] ■ Within the optimal choice of 
parameters, \a\ 2 ~ sinh 2 r ~ n/2 ^s> 1, we fix a total 
number of particles, Nt — pn, distributed in ensembles 
of p independent measurements. There are two concur- 
ring behaviors contributing, in average, to A9: for small 
p (large n) we are in a pre- asymptotic regime charac- 
terized by large oscillations of Y^=i + Nj \ which 
still provides sub shot-noise but not the Heisenberg limit. 
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FIG. 3: (Color online) Demonstration of the Heisenberg limit 
AO ~ 1/Nt- A) Circles are the phase sensitivity obtained 
with the Bayesian analysis as a function of the number of 
measurements p with fixed total number of particles Nt = pn. 
The optimal value, p opt = 30, corresponds to the minimum 
of AO and does not depend on n. Dashed lines are guides to 
the eye. B) Corresponding optimal sensitivity as a function of 
Nt- The dashed line is the asymptotic limit AO — 7. 12 /Nt, 
the solid line is A0 ~ . Shot-noise has been 

y/Popty |a] 2 +sinh 2 r 

included for comparison (dot-dashed line). 



For larger values of p, we saturate the Fisher information 
and obtain A6 = ^Jp/Nx- The prefactor ^fp arises from 
the statistics of independent measurements. As shown in 
figure (|3la), the optimal value is p opt ~ 30. The crucial 
point to notice is that p opt does not depend on Nt- If 
it would, we could not claim the Heisenberg limit. The 
phase sensitivity calculated at p op t is plotted in figure 
(Ob) as a function of Nt (circles). The dashed line is 
the Heisenberg limit A8 = 7.12/Nt, while the solid line 
is Ad w r- 1 For comparison, we include 

\ZPopty I a | 2 +sinh 2 r 

in the figure the shot-noise limit (dot-dashed line). 

We emphasize that an enhancement of phase sensi- 
tivity can be obtained also when only one output port 
is monitored (reduced MZ configuration). A numerical 
calculation of the Fisher information for \a\ 2 ~ sinh 2 r 
shows a strong dependence on 9, the optimal working 
point being close to or tt, depending on the port which 
is monitored. Even if we were not able to numerically 
investigate large values of n, we have strong evidences 
that, asymptotically in n, we obtain a phase sensitivity 
A9 ~ 1/Nt, with a prefactor larger than the one ob- 
tained with the Mach-Zehnder interferometer. 

Discussion. What is the physics underlying the in- 
crease in phase sensitivity using squeezed vacuum light? 
In [l[ Caves associated sub shot-noise to quadrature 
squeezing. Indeed, under the conditions 9 = w/2, and 
\a\ 2 > sinh 2 r, Eq.© reduces to AO = with the 

quadrature X\ = (b* + b)/2. With squeezed- vacuum 
light AX i = e~ r and we recover Eq.(|2|. Conversely, we 
can understand the saturation at the Heisenberg limit by 
quantum interference effects created by the beam split- 
ter. The key point is to notice that the input squeezed 
state has the components S m = when m is odd. For 
the sake of simplicity, we discuss this problem by fixing 
(post-selecting) a total number of particles N — n. The 
input \il> N ) = Ej^-w/a M N I 2 - v)a\N/2 + fi) b is char- 
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FIG. 4: (Color online) Relative number of particles distri- 
bution P(/n) for A) the input state IV'Jv) with post-selected 
N — n and optimal conditions \a\ 2 — sinh 2 r and B) state 
after the first beam splitter, \^p^ s ). C) Phase distribution 
P{4») obtained after a projection of \iPn' S ) over phase states. 
D) Quantity Pnoon as a function of \a\ 2 /n. The maximum 
is reached at \a\ 2 ~ n/2. Here n — 20. 



acterized by a relative number of particles distribution 
P(fjb) = |A M | 2 , where = for odd values of N/2 - ji, 
see Fig.Qa). This creates a relative number of particles 
distribution after the first beam splitter characterized by 
a mean-square fluctuation of the order of N. In partic- 
ular, the distribution has the largest peaks centered at 
/i = ±N/2, see Fig.flUb), which indicates that the corre- 
sponding quantum state after the beam splitter contains 
a large "NOON" component \NOON) - (| N, 0)+|0, N)). 
Such a distribution is typical of states attaining the 
Heisenberg limit A9 ~ \/N . Intuitively, the phase distri- 
bution, P{4>) obtained by projecting a state with heavily 
weighted components at /j, = ±N/2 over phase states 

10 = E^-A72 e ~"W 2 - ")a\N/2 + u) b , is charac- 
terized by oscillations of frequency 2-k/N. This typ- 
ical structure is illustrated in Fig. ((He) where we plot 
the phase distribution obtained by projecting \tp^ s ) = 
e~ t ~2 J °'\ipif) over \<j)). Finally, it is interesting to no- 
tice that the highest "NOON" component is obtained 
when a 2 = n/2, which precisely corresponds to the op- 
timal conditions discussed in Eq.Q. This is illustrated 
in Fig.Qd), where P NO on = \(NOON\ e -^ J *\ip N )\ 2 is 
shown as a function of \a\ 2 /n. 

Contusions. The discovery that interferometric mea- 
surements can be dramatically improved by non-classical 
light has been crucial for the development of modern 
quantum optics [1, [23| • Several states and strategy have 
been proposed in the literature to beat the shot-noise 
limit. Here we have shown that the oldest of these 
proposals, a linear lossless Mach-Zehnder interferometer 
fed by a coherent<S>squeezed- vacuum light [l[ , can indeed 
reach the Heisenberg limit Eq. ((4|) , but only if the whole 
information included in the measurement of the number 
of particles at the output ports is taken into account. 
This requires a feasible analysis of the interferometric 
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data which is provided, for instance, by a Bayesian pro- 
tocol. Moreover, we have also shown that the phase sen- 
sitivity is independent from the true value of the phase 
shift for arbitrary values of squeezing. 
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